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DESCRIZIONE 

a corredo di una domanda di brevetto per invenzione industriale avente per 
titolo: "Metodo per il processamento di immagini con colture di neuroni e 
dispositivo per la sua attuazione" 

a nome: S.I.S.S..A. Scuola Internazionale Superiore di Studi Avanzati 

inventori: Maria Elisabetta Ruaro, Paolo Bonifazi e Vincent Torre 

* * * 

La presente invenzione concerne un metodo per il processamento di immagini 
con colture di neuroni e un dispositivo per la sua attuazione. 
I dispositivi standard al silicio per computer stanno raggiungendo misure 
sempre piCi piccole. Per incrementare il potere di elaborazione dei computer 
attuali, e necessario sviluppare nuovi dispositivi. Sono gte stati costruiti 
dispositivi che contengono molecole organiche 1 " 6 di dimension! inferiori al 
nanometro. Essi risolvono, molto efficacemente, problemi seriali ma, 
nonostante la loro considerevole velocity, sono meno indicati per risolvere 
problemi di parallelismo massivo, come quelli inerenti all'intelligenza artificiale, 
alia visione artificale ed alia robotica ™.Tali dispositivi sono meno indicati per 
parallelism! massivi, a causa della difficolte a formare le numerose 
connessioni tra le varie units di calcolo, per le quali invece i neuroni biologici 
sono ideali. II cervello umano risolve questi problemi, usando element di 
calcolo lenti e spesso inaffidabili 9 " 11 , owero i neuroni, che pero operano 
naturalmente in parallel©. Sono anche stati sviluppati strumenti ibridi, composti 
da neuroni biologici ed elementi di metallo e/o silicio 1SM9 f con lo scopo di 
ottenere neurocomputers, mezzi ideali per risolvere i problemi di parallelismo 
massivo. 



Gli autori della presente invenzione hanno sorprendentemente trovato che, 
attraverso la mappatura di immagini in colture di neuroni cresciute su "array" o 
matrici di multielettrodi disponibili commercialmente (MEA), e' possibile 
elaborare in pochi second! immagini digitali, indipendentemente dalla loro 
dimensione, contrariamente a quanto awiene nei computer seriali 
convenzionali. L'invenzione mostra la fattibilite e I'utilite dei neurocomputer e 
stabilisce la base per una nuova generazione di strumenti per computer. 

Forma pertanto oggetto della presente invenzione un dispositivo per il 
processamento in parallelo di immagini comprendente almeno una "multi 
electrode array" MEA, su cui 6 cresciuta una coltura omogenea di cellule 
neuronal!. La MEA viene utilizzata sia per la stimolazione dei neuroni sia per 
la registrazione della loro attivita elettrica. 

E' oggetto dell'invenzione un metodo per il processamento in parallelo di un 
immagine digitale comprendente le seguenti fasi: 

a) mappatura di una immagine digitale (li, 2 (x,y)) (INPUT) con 1 o 2 bit di 
risoluzione (li(x,y)) o l 2 (x,y)) a seconda che si abbia un'immagine 
rispettivamente di 1 o 2 bit) di NxN pixel in impulsi di voltaggio di 2 o 4 livelli di 
intensity applicati ad una matrice NxN di elettrodi integrati su un piattino 
(MEA) su cui sono mantenute in coltura cellule neuronali connesse 
spontaneamente fra loro a formare una rete; 

b) elaborazione deirimmagine da parte della coltura neuronale con il nucleo o 
kernel di convoluzione: 

h(p,a,t) = A(t) exp((p - p(t)) / 2 a(t) 2 ) (1 ) 



Questa elaborazione d compiuta dalle connessioni che si stabiliscono trai 
neuroni dissociati a seguito delta messa in coltura dei neuroni. Questa e una 
componente essenziale del nuovo dispositivo poiche ne costituisce il 
meccanismo di calcolo; 

c) registrazione dell'attivita elettrica della coltura neuronale a mezzo di 
registrazione con la matrice di multielettrodi dei segnali elettrici (in voltaggio) 
extracellulari e rivelazione, per ogni singolo elettrodo e in successivi intervalli 
di tempo, dei segnali o span associati alia generazione dei potenziali di azione 
dei neuroni in coltura. 

Preferibilmente si registra la frequenza di sparo FR(x,y,t) (OUTPUT) misurata 
dall'elettrodo in posizione (x,y) ed in un intervallo di tempo centrato in t. 

Preferibilmente, la relazione tra INPUT ed OUTPUT sono legati fra loro 
daU'equazione: 

FR(x,y,t) = k 2 (x,y)** hfoctf) (2) 
dove ** indica una convoluzione bidimensionale. 

U equazione (2) descrive come e' processata una singola immagine digitale di 
1 o 2 bit (l 1(2 ). 

Nel caso in cui Timmagine digitale in INPUT sia definita con 8 bit (Ufay)), essa 
viene scomposta altemativamente in 4 o 8 immagini rispettivamente da 2 o 
1bit (l mi ) dove m 6 uguale rispettivamente a 2 o 1 secondo I'equazione: 

8/m 

Wx^LU^ (3) 
1 

ed ogni singola immagine l m j viene processata indipendentemente con il 
metodo sopra descrltto. Queste vengono poi ricomposte in una unica 



immagine ad 8 bit e I'intero processo (scomposizione, filtraggio, 
ricomposizione) e' descritto dall'equazione: 

8/m 

Z|U2 m(M) ** h (p ja|t ) (4) 
1 

il cui risultato finale 6 rimmagine l B (x t y) processata con una risoluzione ad 8 
bit. 

L'invenzione verra ora descritta secondo forme esplicative ma non limitative, 
facendo riferimento alle seguenti figure: 

Figura 1: Mappaggio di un'immagine nella stimolazione di una coltura 
neuronale. a) immagine digitale binaria 6x10 di una L usata come indice di 
stimolazione di una coltura neuronale cresciuta su una MEA 6x10 prodotta da 
MCS. b) coltura neuronale ottenuta dalla dissociazione di neuroni da 
ippocampo (vedi Metodi). Un ingrandimento della coltura di neuroni sull'area 
della MEA evidenziata dalle lettere B e C (rettangolo bianco) e mostrata 
nell'inserto. c) Attivita elettrica neurale evocata e registrata dalla MEA. 
L'elettrodo indicato dalla freccia § stato usato come terra, d) tre registrazioni 
rappresentative di attivita'evocata da stimoli di voltaggio di 0.3, 0.6 e 0.9 V. e) 
AFR (vedi metodi) registrati da un elettrodo rappresentativo (lo stesso usato in 
d) a stimolazioni differenti di voltaggio, come indicato nel pannello. f) AFR con 
stimolazioni ripetute a intervalli differenti come indicato nel pannello. Dati 
ottenuti da 50 diversi esperimenti della stessa stimolazione. Tempo 0 
corrisponde airapplicazione dell stimolazione. In e e f i conteggi sono stati 
effettuati su intervalli di 10 millisecondi. 

Figura 2. Diffusione dell'eccitazione attraverso la coltura neuronale. a, b,c,d 
mostrano la AFR negli intervalli di tempo di 1-6, 4-9, 7-12 e 12-17 msec 
rispettivamente dopo la stimolazione della fila superiors deglt elettrodi con un 



6 

impulso di voltaggio di 0.6 V. I punti rossi sono i dati sperimentali e linee solide 
sono le curve teoriche ottenute con I'eq (1), In a e b la curva e stata ottenuta 
fissando paOea rispettivamente a 890 e 1240. In c e d, le curve sono state 
ottenute fissando p a 920 e 1750 e a a 980 e 1130 rispettivamente, e-h 
immagini processate dalla coltura neuronale nei rispettivi intervalli di tempo H: 
filtraggio digitale di un'immagine binaria, da 6x10 pixel, di una barra 
orizzontale (riga in alto) nei rispettivi intervalli di tempo. II filtraggio e stato 
ottenuto con i valori delle curve teoriche riportati in a,b,c,d. AFR in e-l sono 
rappresentati secondo la mappa di colore (vedi sezione Metodi) riprodotta 
nella figura. 

Figura 3. reproducibility del filtraggio deil'immagine da parte delta coltura 
neuronale. A: Attivite elettrica neuronale extracellulare registrata da uno 
stesso elettrodo in risposta a tre identiche stimolazioni. b,c,d f e,f: AFR e CV 
ottenuti su una finestra di 10 msec corrispondenti agli elettrodi posti a distanza 
crescente rspetto alia barra di stimolazione e indicati con la stessa 1 lettera 
mostrata neirinserto. g: ciascuna fila riproduce immagini ottenute da una 
singolastimolazione, nei quattro intervalli di tempo, indicati sopra ciascuna 
colonna. La FR in ciascuna immagine e' rappresentata secondo la mappa di 
colon (vedi sezione Metodi)jiportata a destra della figura. 
Figura 4: Filtraggio passa banda della coltura neuronale. a,b: Profilo e sua 
trasformata di Fourier in risposta alia stimolazione neirintervallo di tempo tra 
1 e 6 msec. I punti rossi sono i punti sperimentali riportati in Fig. 2. La linea 
continua da eq (1) con oep uguali a 890 e 0 rispettivamente. c,d: come in a,b 
ma neirintervallo di tempo tra 4 e 9 msec. La curva in c 6 ottenuta dall'eq (1) 
con a e p uguali a 1240 e 0 rispettivamente. e.f: risposta all'impulso e sua 



trasformata di Fourier ottenute come differenza traaecetrabed 
rispettivamente. In a f c,e nelle ascisse 6 usata la stessa unite arbitraria (a.u.). 
Similmente in b,d,f. In b,d,ed e, F1(k) e F2(k) sono la trasformata di Fourier 
della cx>rrispondente Gaussiana in a,c; con k uguale a 2tcA,, dove X 6 la 
lunghezza dell'onda spaziale della trasformazione Fourier. g,h: pannelli di 
sinistra: filtraggio passa banda della coltura neuronale delle due immagini 
binarie, composte da una barra orizzzontale (-) e una L, rispettivamente, 
ottenute sottraendo la AFR nell'intervallo di tempo di 1-6 e 5-10 msec; pannelli 
di destra: filtraggio digitate ottenuto convolvendo Timmagine binaria originale 
con la differenza delle due Gaussiane mostrate nel pannello e. Le barre sottili 
indicano gli elettrodi stimolati. La codificazione del colore come descritta nella 
sezione Metodi e riprodotta alia destra della figura. i,l: filtraggio passa basso di 
due different! immagini a 8 bits. Pannelli di sinistra: immagini originali a 8 bits. 
Pannelli centrali: filtraggio passa basso delle immagini, ottenuto con le colture 
neuronali. Pannelli di destra: filtraggio passa basso delle immagini, ottenuto 
da una convoluzione digitale dell'immagine originale a 8 bits, con il profilo 
Gaussiano mostrato in 4a. La codifica del colore come descritta nella sezione 
Metodi 6 riprodotta a destra della figura. 
METODI 

Preparazione della coltura dei neuroni 

I neuroni primari da ippocampo di ratti Wistar di 3giorni, sono ottenuti 
essenzialmente come descritto 32 . La dissezione viene fatta in 200 pM acido 
kynurenico e 25 pM 2-amino-5-fosfonovalerato (APV) (Sigma, St. Louis, MO) 
in ghiaccio. Le cellule vengono dissociate meccanicamente facendo passare I 
frammenti di ippocampo per 20-30 volte attraverso una punta Gilson blu. II 
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terreno di coltura e' il minimal essential medium (MEM) con sali di Earie e 
Glutamax I (Life Technologies, Gaithersburg, MD), con aggiunta di 5% di 
siero bovino fetale, 7 mg/ml D-glucosio, 3.6 mg/ml HEPES, 0.1 pg/ml biotina, 
1.5 pg/ml vitamina B12, 30 pg/ml insulina.e 100 pg/ml transferrina bovina. 
II rivestimento della superficie della MEA con molecole di adesione viene fatto 
incubando per una notte a 37°C con 1 ml di soluzione di poliornitina alia 
concentrazione di 50//g/ml di acqua. Le piastre sono successivamente 
asciugate alFaria. 20 minuti prima della piastratura la superficie della MEA 
viene ulteriormente coperta con uno strato sottile di Matrigel (Beckton- 
Dickinson) solo sulla regione della matrice dl elettrodi. 

100 fj\ di sospensione cellulare alia concentrazione di 8x1 0 5 cells/cm 2 . 6 
depositata sulla zona degli elettrodi di un MEA prerivestito con fattori di 
adesione Le cellule sono lasciate a temperatura ambiente per 20 minuti per 
permetterne la deposizione sul fondo della MEA, quindi viene aggiunto 1 ml 
del terreno di coltura e la MEA viene incubata in un'atmosfera di 5% C0 2 a 
37°C. La proliferazione di cellule non neurali viene bloccata con Taggiunta 
di2.5-5.0 pM citosina-D-arabinofuranoside dal secondo giomo di coltura in poi. 
Meta del terreno 6 cambiato due volte alia settimana. Le registrazioni sono 
effettuate dalla 3a settimana dopo il piastramento fino a 3 mesi dopo. La 
stessa piastra pud essere utilizzata per le registrazioni elettriche piu volte. 
Durante le registrazioni elettriche, le piastre sono coperte da uno speciale 
tappo ottenuto da MultiChannelSystem (MCS), che riduce gli scambi gassosi 
e impedisce I'evaporazione e la contaminazione. 
Registrazioni elettriche e stimolazione dell'elettrodo 
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II slstema commerciale fornito dalla MCS (www.multichannelsystems.com) § 
usato per la stimlolazione e la registrazione elettrica. I risultati riportati si 
riferiscono ad un matrice di microelettrodi 6x10, con gli etettrodi distanziati tra 
loro di 500 micron. Questa matrice viene usata per la stimolazione e la 
registrazione elettrica. 

Attraverso contatti in oro f la matrice di elettrodi § connessa ad 

preamplificatore/filtro/amplificatore a 60 canali, la cui ampiezza di banda e' 

i 

10Hz-3kHz. (MEA 1060-AMP) che riorienta i segnali verso un ulteriore 
processo elettronico (ad es. ampliflcazione e conversione AD). Le acquisizioni 
del segnale sono gestite da un controllo software. Un termostato (HC-X), 
posto al di sotto del MEA mantiene la temperatura a 37°C, . II MEA fomito da 
MCS 6 in grado di digitare in tempo reale a 20 kHz tutte le registrazioni di 
voltaggio Vij ottenute dai 60 elettrodi metallici. Un elettrodo e usato come terra 
(Fig. 1c). I dati campionati sono trasferiti in tempo reale ad un disco fisso per 
un processamento successivo. Ciascun elettrodo metaliico pud essere usato 
per la registrazione o per la stimolazione in quanta il sistema MCS non 
consente un passaggio controllato da computer da stimolazione a 
registrazione . Le stimolazioni di voltaggio Sij consistono in impulsi bipolari di 
durata di 100 msec a ciascuna polarity, per un'ampiezza variante da 200 mV 
a 1 V, iniettata attraverso il generatore di stimoli STG1004 . Un artefatto della 
durata 5-20 msec causato dalla stimolazione elettrica e indotto sugli elettrodi 
di registrazione. Questo artefatto e rimosso dalle registrazioni elettriche 
durante le analisi dei dati. 
Analisi dei dati 
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I dati acquisiti sono stati analizzatl con il software MatLab (The Mathworks, 
inc.). 

Rimozione dell'artefatto: L'artefatto a ciascun elettrodo e per ciascun schema 
di stimolazione e stimato e sottratto dalle registrazionl del voltaggio. L'artefatto 
e stimato nel seguente modo: per ciascun schema di stimolazione e a ciascun 
elettrodo la risposta del voltaggio come media di tutte le prove (tipicamente 
50) e calcolata e interpolata da 2 polinomi di 9 grado. I due polinomi 
dispongono i dati in un intervallo di tempo di 0.5-25 ms e 7.5-100 ms dopo la 
stimolazione. II primo polinomio e usato per valutare l'artefatto nell'intervallo di 
tempo da 0.5 a 7.5 msec, mentre il secondo neH'intervallo di tempo da 7.5 e 
82.5 msec. L'artefatto, cosl valutato, e sottratto dal segnale di voltaggio 
originate. 

Calcolo della frequenza di sparo (firing rate) (FR): Vij e il voltaggio registrato 
all'elettrodo (i,j) e oij e la deviazione standard del rumore calcolato in un 
periodo di almeno 1 sec in cui nessun picco o sparo di neurone e visivamente 
osservato. Questo indice di frequenza di sparo FRij(t) al tempo t=(t1+t2)/2 e il 
numero di tutti gli incroci dei livelli Vij al di sopra di una soglia fissata a 5Vij 
computata in un intervallo di tempo tra t1 e t2. Questo FRij(t) conta gli spari 
prowenienti da differenti neuroni, che fanno un buon contatto elettrico con 
I'elettrodo (ij). La a ij del rumore varia per ciascun elettrodo da 3 a 6 //V. 
L'indice di frequenza di sparo medio AFRij(t) e calcolato dalla media di FRij(t) 
di tutte le prove. AFR(t) e calcolato su un intervallo di 10 msec. 
Elaborazi one dell'immaqine 

I MEA con almeno piu di 56 elettrodi che fomiscono registrazioni elettriche di 
picchi o spari chiari sono usate per I'elaborazione dell'immagine. Data 
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un'immagine /ij di N x N pixels ed un MEA con N x N elettrodi, il livello di grigio 
di un pixel (i j) di / e convertito in un appropriata stimolazione di voitaggio Sij 
deH'elettrodo (ij). II MEA fornisce i segnali di voitaggio Vij composti da 
potenziali di azione o spari prodotti dai neuroni. II processamento 
deirimmagine /ij 6 il set di risultati FRij (t), cosicch6 a tempi differenti t, c'§ un 
differente processamento dell'immagine originate / ij. 

Mappaggio lij in Sij. Sia V1/2 la stimolazione di voitaggio che evoca la mete 
della AFR massimale 10 msec dopo I'impulso di voitaggio. Se / ij £ 
un'immagine binaria, ad es. se i suoi livelli di grigio sono sia 0 che 1 , allora Sij 
sara 3/2 * V1/2 se / ij e 1, 0 nell'altro caso. Se / if e un'immagine a 2 bits, ad 
es. se i suoi livelli di grigio sono 0, 1 ,2 o 3, allora Sij sara 0 se / ij 6 0, Sij sara 
1/2*V1/2se//y e1, Sij saraV1/2 se//y e2e Sij sara3/2*V1/2 se I ij e 
3. 

Riempimento di elettrodi silenti e lisciamento. Quando un elettrodo (ij) e 
silente, cioe quando non si registra alcun picco o sparo, il buco corrispondente 
nelPimmagine processata e riempito assegnando a FRij(t) il valore ottenuto 
calcolando la frequenza di sparo medio degli elettrodi confinanti - ad es. 
elettrodi a una distanza di of 500 //m. FRij(t) degli elettrodi stimolati e 
determinato attraverso I'estrapolazione dagli elettrodi attivi confinanti usando 
la eq(1). Tutte le immagini che sono state processate ed elaborate dalla 
coltura neuronale avevano al massimo 3 elettrodi silenti, incluso quello usato 
come terra. Di conseguenza si usavano solo colture su MEA in cui quasi tutti 
gli elettrodi erano funzionanti. Spesso, il valore di FRij(t) e stato lisciato 
arrotondato il suo valore facendo la media sugli elettrodi confinanti (i-1 j), 
0+1 J). (ij-1)e(ij+1). 
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Processamento di immagini a 8 bits. Le immagini a 8 bits / 8 (x,y) sono 

scomposte come: 

H ln(x.y)2™ 
1 

dove In(x,y) 6 un'immagine a 1 bit, o 

4 

* Ii2(x,y) 2 2(M> 
1 

dove l, 2 (x,y) S un'immagine a 2 bits. Le 8 immagini l h (x f y) a 1 bit o le 4 
immagini li 2 (x,y) a 2 bit sono processate come descritte sopra e il loro risultato 
§ sommato come descritto in eq (4). 

Risultato delta codifica del colore. Le immagini processate FRij(t), AFRij(t), o 
la loro combinazione (per filtraggio passa basso e/o processamenti a 8 bits) 
sono mostrate usando una procedura standard di codifica del colore: per 
immagini filtrate passa basso, i valori di FRij(t) o AFRij(t) sono normalizzati ad 
avere valori tra 0 e 1 dividendo per il valore corrispondente massimale tra 
tutti gli elettrodi nell'intervalio di tempo 0-25 ms e la mappa di colore (di 256 
colorij e normalizzata tra -1 e 1. Per immagini filtrate passa banda ottenute 
come differenza di FRij(t) o AFRij(t) a due tempi differenti - il risultato 
derivante 6 normalizzato tra -1 e +1, dividendo per il suo~valore assoluto 
massimo e la corrispondente mappa di colore (di 256 colori) e normalizzata tra 
-1 e 1. La codifica del colore di immagini filtrate digitalmente e ottenuta con la 
normalizzazione dei valori di intensity tra -1 e 1 ed una simile codifica del 
colore. Pertanto, con questa codifica, il processamento delle immagini a 1,2 o 
8-bits ha la medesima mappa, cio§ il risultato ha 256 colori differenti. Se 
desiderata si pud utilizzare una diversa codifica. 
RISULTATI 
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I MEA con 60 o piCi elettrodi possono essere ottenuti da centri di ricerca o 
acquistati cxxnmercialmente. I MEA sono fabbricati con diverse geometrie di 
elettrodi, tra cui una griglia quadrata regolare con spazi tra gli elettrodi che 
variano da 50 a 500 //m. Gli elettrodi individual! sono solitamente coperti da un 
sottile strato di platino e hanno forma quadrata con lati da 10 a 30 jjm. La 
grande maggioranza dei MEA disponibiii attualmente e le matrici di pixels per 
camere CCD condividono la stessa geometria di una griglia quadrata. Questa 
osservazione ha ispirato il disegno di un dispositivo per la elaborazione di 
immagini dove I'elemento di computo 6 una coltura di neuroni cresciuto su 
una base MEA: I'immagine e' mappata alio stimolo voltaico di una cultura di 
neuroni e I'attivita elettrica evocata e* presa come segnale di uscita del nuovo 
dispositivo. 

Descrizione dettaaliata del dispositivo. 

Un'immagine di /ij of N x N pixels (Fig. 1a) con Tusuale geometria quadrata e 
codificata nel segnale di entrata di un MEA con N x N elettrodi (Fig. 1b), cosi 
che il livello di grigio di pixel (ij) di /ij e convertito in una appropriata 
stimolazione di voltaggio Sij delPelettrodo (ij) (si vedano i Metodi per una 
descrizione dettagliata di questa codifica). II segnale di uscita del dispositivo e 
composto dai segnali di voltaggio Vij registrati con il MEA (Fig. 1c). Questi 
segnali sono composti da potenziali di azione o spari prodotti dai neuroni nella 
coltura (mostrato nell'inserto del pannello b), e le loro statistiche sono usate 
per ottenere una codifica dell'immagine processata. PiCi specificatamente, il 
segnale di uscita del dispositivo e FRij(t), cioe la frequenza di sparo di tutti i 
neuroni registrati all'elettrodo (i j) nel intervallo di tempo tra t - At et + At, 
successivo alia stimolazione dell'elettrodo al tempo 0. In questo modo, per 
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ciascuna stimolazione Sij codificante per I'immagine /ij, si ottiene una serie di 
segnali di uscita FRij (t) rappresentanti r elaborazione di /ij al tempo t. 
In quello che segue sara mostrato che una coltura di neuroni da ippocampo di 
ratto (vedi Metodi) pud essere usata per ottenere un filtraggio passa basso e 
passa banda di immagini digitali. 

L'intervallo dinamico, il tempo di ciclo e la reproducibility della risposta del 
dispositivo proposto sono stati analizzati stimolando una fila di elettrodi del 
MEA con lo stesso impulso di voltaggio. Questa stimolazione corrisponde ad 
un'immagine nera con una barra sottile chiara. Sono stati usati brevi voltaggi 
bipolari che durano 200 psec e la loro ampiezza e'stata progressivamente 
aumentata da 300 mV to 900 mV. I picchi registrati da un elettrodo 
aumentavano in frequenza e spesso apparivano picchi di forma diversa, 
prodotti da un differente neurone (Fig. 1d). La frequenza di sparo media (AFR) 
di tutti i picchi rilevati (vedi Metodi) aumenta con la stimolazione del voltaggio 
(Fig. 1e). L'intervallo dinamico della AFR, tuttavia, era piuttosto stretto: 
solitamente nessun picco era evocato con impulsi di voltaggio al di sotto di 
200 mV e, una risposta massimale saturante era evocata con un voltaggio di 
stimolazione di circa 1 V. Nella grande maggioranza degli esperimenti e stato 
possibile distinguere in modo affidabile quattro livelli di attivite evocata, 
indicanti che la coltura neuronale pud codificare 2 bit 

Per determinare il tempo del ciclo del dispositivo, la stessa stimolazione e 
stata ripetuta ad intervalli da 100 msec a 10 secondi. Usando un intervalio 
maggiore di 1 o 2 secondi tra stimolazioni successive , la AFR aveva due 
component!: uno, evocato con un ritardo di pochissimi msec e di circa 15 
msec di durata, seguito da un secondo componente di durata di circa 100 
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msec. L'ampiezza del primo componente non era signlficativamente 
influenzata dalla diminuzione del tempo di ripetizione da 10 sec a 0.1 msec 
(Fig. 1f). L'ampiezza del secondo componente era significativamente diminuita 
a brevi tempi di ripetizione ed era stabile per tempi di ripetizione superiori ai 4 
secondi. Pertanto, il nuovo dispositivo pu6 processare 2 bits con un tempo di 
ciclo che varia da 0.25 a 10 Hz, a seconda se 6 considerate il primo o il 
secondo componente nella risposta. 

La coltura neuronale cresciuta sul MEA costituisce un network bidimensionale 
e le sue propriety di filtraggio sono meglio analizzate utilizzando una lunga 
barra come stimolo spaziale. In questo modo, data una coltura omogenea, la 
caratterizzazione di un network bidimensionale e ridotta alia comprensione di 
un problema monodimensionale molto piu semplice: i sei elettrodi della fila 
superiore sono stati usati per la stimolazione e la AFR evocata nei nove 
elettrodi della stessa colonna e stata misurata. 

Nei primi momenti, come ad esempio neirintervallo di tempo tra 1 e 6 msec 
(Fig 2a) I'AFR evocato decade spazialmente attraversando la coltura 
neuronale come una funzione Gaussiana con una deviazione standard a di 
890 /vm, corrispondente a 1 .8 pixels. Nell'intervallo di tempo tra 4 e 9 msec 
(Fig 2b) Pattivita' elettrica decade con una funzione Gaussiana similare ma 
con una a t corrispondente a1240//m corrispondente a 2.5 pixels. Dopo circa 
10 msecs il picco dell'AFR evocata si allontana dagli elettrodi stimolati ed e' 
descritta come una funzione Gaussiana centrata ad una distanza p dagli 
elettrodi stimolati corrispondente a 2.5 pixels. Dopo ulteriori 10 msecs il picco 
deTAFR evocata si muove dagli elettrodi stimolati (Fig 2d). II processamento 
della coltura neuronale dell'immagine /ij composta da una barra chiara di 6 
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pixels nella parte superiore e rappresentata dall* immagine codificante la FR, 
come mostrato nei pannelli e, f. g, h. Nei pannelli i e j sono mostrate le 
convoluzioni digitali dell" immagine iniziale /ij con un kernel Gaussiano con 
una deviazione standard a di 1.8 e di 2.5 pixels rispettivamente._Nei pannelli i 
e j sono mostrate le convoluzioni a tempi successivi di 9.5 e 14.5 msec. FRij 
(t) nei pannelli g e h rispettivamente e uno spostamento del filtraggio passa 
basso deH'immagine originale. 

Gli esperimenti dove una fila (o una colonna) di elettrodi, o elettrodi individuali, 
sono stimolati indicano che il processamento spaziale-temporale della coltura 
neuronale e descritto da una risposta di impulso radiale. 

h(p,a,t) = A(t) exp«p - p(t)) / 2 a(t) 2 ) p^+y 2 (1 ) 

cioe una usuale funzione Gaussiana usuale , centrata su p(t) e con una 

varianza a 2 (t) che varia nei tempo . Pertanto, data un immagine a 2 bit / 2 bit 

(x,y), la risposta di uscita del dispositivo varia con il tempo ed e: 

FR (x,y,t) =/ 2 bit (x.y)** h(p,a,t) 

in cui p e a dipendono dal tempo e ** indica una convoluzione bidimensionale. 
L'immagine in Fig.2k, ottenuta per convoluzione di /ij con il kernel (1) avente i 
valori di 1250 pm e 850 //m per p e a rispettivamente, mostra le stesse 
caratteristiche di quelle in Fig.2g che rappresenta il filtraggio neuronale 
osservato sperimentalmente. 

Diversamente dagli strumenti in silicio, i neuroni biologici sono affetti da un 
significativo rumore ed e pertanto essenziale valutare I'affidabilita del 
dispositivo proposto. La Fig.3a mostra tre differenti risposte elettriche alia 
stessa stimulazione: il numero del picchi o span evocati era leggermente 
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variabile, ma i primi picchi o span evocati awenivano con una variabilita 
molto piccola di circa 1 msec. 

Al picco della risposta evocata, il coefficiente di variazione (CV) delia 
frequenza di sparo p, calcolato su una finestra di 10 msec - era sempre 
piccolo, circa 0.2. Alia fine il CV aumentava fino a circa 1 e anche piCi. Con 
finestre piu' larghedi 50 msec, il CV era sempre inferiore a 0.5. II CV della 
risposta evocata non si modificava in modo significativo con la distanza 
dell'elettrodo di registrazione dal sito della stimolazione (Fig. 3c,d,e,f). La 
Figura 3g illustra le immagini ottenute in singoli esperimenti o prove - senza 
essere mediati - quando la fila superiore degli elettrodi e stimolata con lo 
stesso impulso di voltaggio di 600 mV. All'inizio (prima e seconda colonna) le 
immagini processate erano abbastanza simili. Dopo i primi 15 msec, le 
immagini processate differiscono da prova a prova, in accordo con Palto CV 
delle registrazioni elettriche. Questi dati mostrano che immediatamente dopo 
la stimolazione del voltaggio c'e una "buona" finestra temporale durante la 
quale iP elaborazione della coltura neuronale e riproducibile, dando luogo ad 
un computo affidabile. 

Dato il filtraggio spazio-temporale di eq (1), durante i primi millisecondi 
successivi alia stimolazione elettrica e possibile attuare molto velocemente un 
filtraggio passa basso e passa banda di un'immagine. NeH'intervallo di tempo 
di 1-6 msec, la risposta all'impulso della coltura neuronale 6 una funzione 
Gaussiana con un a di circa 890 /vm (Fig. 4a), ma 2 o 3 msec pio; tardi il 
valore di a aumenta a 1240 //m (vedi fig 4c), i due filtraggi sono passa 
basso (vedi la loro densitO di frequenza di energia nei pannelli 4b e 4d), ma la 
loro differenza (Fig. 4e) e' a passa banda (Fig 4f). II filtro passa banda delle 
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immagini binarie rappresentanti caratteri semplici come un - e L ottenuta con 
la coltura neuronale sono mostrate nelle Fig. 4g,h. II confronto con lo stessb 
flltro effettuata da un computer digitale mostra che il dispositivo opera come 
desiderata. 

La coltura neuronale puo' essere utilizzata anche per elaborare immagini 
digital! a 8 bits. Si consideri l 8 (x,y) come un'immagine con livelli grigi di 8 bits 
nella posizione (x,y). Allora l 8 (x,y) pud essere scomposta come: 
a/m 

l 8 (x.y)=Zi l m ,(x,y)2 m(M > (3) 
1 

dove l mt (x,y) e una quantita con 1 o 2 bits (con m uguale rispettivamente 
quindi a 1 o 2). Data questa scomposizione, la elaborazione di un'immagine a 
8 bits 6 ottenuto come: 

8/m 

I, U 2 m(M > ** h(p,o,t) (4) 

Elaborando con la cultura neuronale, indipendentemente, le 4 immagini a 2 
bits o ie 8 immagini a 1 bit, si ottiene un filtraggio passa basso o passa banda 
di un'immagine a 8 bits, come quella mostrata nella Fig. 4i. I, pannelli di 
sinistra. Un filtraggio passa basso delle immagini originali a 8 bit, ottenute 
dalla coltura neuronale e dal filtraggio digitale con una funzione Gaussiana, 
sono mostrate nei pannelli di centro e di destra, rispettivamente. La alta 
similarity dell'immagine nel pannello di centro e di destra mostra che il 
dispositivo proposto pud processare anche immagini a 8 bits. 
E' pertanto dimostrato che, crescendo colture neuronali su multi electrode 
arrays (MEA), si ottiene un nuovo dispositivo di computo composto da neuroni 
biologici ed elettrodi metallici. I due maggiori viintaggi di questo dispositivo 
sono: prima di tutto, un filtraggio passa basso e passa banda pu6 essere 
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ottenuto in meno di 10 msec, indipendentemente dalla dimensione 
deirimmagine da processare. Come conseguenza il dispositivo e 
potenzialmente in grado di processare ampie immagini digital! (2000 x 2000 
pixels) piti velocemente della magglor parte dei computer digital! noti; in 
secondo luogo, Puso dei neuroni biologici come strumenti di compute apre una 
nuove prospettive nelle quali I'informatica e I'ingegneria dei calcolatori pud 
attingere dall'esperienza e dalla tecnologia della biologia cellulare e 
dell'ingegneria genetica. 

futility del dispositivo proposto dipende dalla dimensione del processamento 
parallelo. II dispositivo proposto non e conveniente per la elaborazione di 
piccole immagini, le quali possono essere accuratamente elaborate da 
computer digitali standard. Diventa utile, fornendo migliori prestazioni dei 
computer digitali, quando devono essere processate immagini grandi con piu' 
di 1000x1000 pixels. II processamento di queste immagini, tuttavia, richiede 
lo sviluppo di MEA con piu di 1 milione di elettrodL Inoltre, lo sviluppo di MEA 
con un numero molto alto di elettrodi e la soluzione di tutti i problemi di 
interfaccia, Tuso efficiente del neurocomputer richiede anche un appropriate 
quadro di compute. I neuroni biologici sono lenti e non molto affidabili come 
element! di compute, ma lavorano naturalmente in parallelo. Sono ideali per la 
soluzione di problemi di parallelismo massivo, quando I'affidabilite di un 
singolo elemento di compute non e critica. I neuroni biologici e probabilmente 
tutti i neurocomputer non sono in grado di imitare una macchina Turing 20 9 cioe 
un dispositivo di compute seriale e preciso. Una delle maggiori attrazioni dei 
neurocomputer e la possibility di usare la intrinseca adattabilite dei neuroni 
biologici, derivati da milioni di anni di evoluzione. La possibility di sfruttare il 
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potenzlamento (LTP)21 e la depressione (LTD)22 a lungo termine nelle 
colture neuronali fomisce una naturale implementazione di algoritmi elaborati 
per reti neurali artificiali (ANN)23-24. E' anche possibile usare la tecnologia 
delle cellule staminali25-27 per ottenere popolazioni di neuroni con proprieta 
speciflche e rllascianti neurotrasmettltorl selezionati cos! da avere colture 
neuronali con rapporti controllati di neuroni inibitori ed eccitatori, possibilmente 
anche con un'affidabilite piCi alta degli elementi di computo individuali. La 
possibility di direzionare la crescita dei neuroni lungo direzioni spaziali 
specifiche28-30 consentira la fabbricazione di un ampia varieta di filtri spaziali, 
imitando le proprieta recettive di campo di neuroni nelle aree visive 31. 

II dispositivo ha una programmabilita limitata. II calcolo rappresentato dal 
dispositivo pud essere controllato da differenti chimici aggiunti al mezzo 
extracellulare ma non e verosimile ottenere la versatility dei soliti computer 
digitali. 

Questo dispositivo e destinato ad essere alia base di una nuova generazione 
di strumentazione di calcolo, sviluppato daila sinergia della scienza dei 
materiali e della biologia cellulare. Questi mezzi di calcolo sembrano avere 
caratteristiche simili a quelle degli esseri umani come I'apprendimento e 
I'adattabillta. 
Biblioarafia 

1. Chen, J., et al. Large On-Off Ratios and Negative Differential Resistance in 
a Molecular Electronic Device. Science 286, 1550-1552 (1999). 

2. Collier C. P., et al. Electronically Configurable Molecular-Based Logic 
Gates. Science 285, 391-394 (1999). 



21 



3. Collier, C. P., et al., A [2]Catenane-Based Solid State Electronically 
Reconfigurable Switch. Science 289, 1 172 (2000). 

4. Lee, I., Lee, J., & Greenbaum, E. Biomolecular Electronics: Vectorial Arrays 
of Photosynthetic Reaction Centers. Phys. Rev. Lett. 79, 3294 (1997). 

5. Reed M. A., et al. Conductance of a Molecular Junction, Science 278, 252- 
254. (1997). 

6. Roth K. M. et al., Molecular approach toward information storage based on 
the redox properties of porphyrins in self-assembled monolayers, J. Vac. Sci. 
Tech.B 18, 2359 (2000). 

7. Marr. Vision: A Computational Investigation into the Human Representation 
and Processing of Visual Information. San Francisco. CA:Freeman & Co., 
(1982) 

8. Rumelhart, D. E. and J. L. McClelland. Explorations in Parallel Distributed 
Processing. MIT Press. Cambridge, Massachusetts. (1988) 

9. Nicholls J.G., et al. From Neuron to Brain: A Cellular and Molecular 
Approach to the Function of the Nervous System 4TH edition. Sinauer 
Associates, Incorporated (2000) 

10. Shadlen MN, Newsome WT (1994) Noise, neural codes and cortical 
organization. Curr Opin Neurobiol 4: 569-579 

11. Zoccolan D, G. Pinato & Torre V. Highly variable spike trains underlie 
reproducible sensory-motor responses in the medicinal leech. J NeuroscL (in 
press ) (2002). 

12. Gross, G.W., et al. A new fixed-array multimicroelectrode system designed 
for longterm monitoring of extracellular single unit neuronal activity in vitro. 
NeuroscL Lett. 6, 101-105(1977) 



22 

13. Pine, J. Recording action potentials from cultured neurons with 
extracellular microcircuit electrodes. J. Neurosci. Methods, 2, 19-31 (1980) 

14. Novak, J.L. & Wheeler, B.C. Recording from tbe aplysia abdominal 
ganglion with a planar microelectrode array. IEEE Trans, Biomed. Eng. 33 196- 
202(1986) 

15. Jimbo Y. & Kawana A. Electrical stimulation and recording from cultured 
neurons using a planar electrode array, Bioelectrochem. Bioenergetics 29, 
193-204(1992) 

16. Martinoia.S., et al. A general purpose system for long-term recording from 
a microelectrode array coupled to excitable cells. J. Neurosci. Methods 48, 
115-121 (1993) 

17. Vassanelli, S., & Fromherz, P., Neurons from Rat Brain coupled to 
Transistors. Applied Physics A 65, 85-88 (1997) 

18. Zeck, G., & Fromherz, P. Noninvasive neuroelectronic interfacing with 
synaptically connected snail neurons immobilized on a semiconductor chip. 
Proc. Nat Acad.Sci. USA 98, 10457-10462 (2001). 

19. Bonifazi,P., & Fromherz, P. Silicon Chip for Electronic Communication 
between Nerve Cells by Noninvasive Interfacing and Analog-Digital 
Processing. Advanced Materials 14, 1 190-1 193 (2002) 

20. Cutland, N. Computability, Cambridge University Press, (1980). 

21. Bliss, T.V.P. & Collingridge, G.L. A synaptic model of memory: long-term 
potentiation in the hippocampus. Nature, 361 31-39 (1993). 

22. Linden, D.J. & Conner, J. A. Long-term synaptic depression . 
Ann.Rev.Neurosci. 18, 319-357 (1995). 



23 



23. Introduction to the Theory of Neural Computation. Hertz J.. Krogh A., 
Palmer R.G., (Santa Fe Institute Studies in the Sciences of Complexity. 
Lecture Notes, Vol 1 ) (1 991 ). 

24. Hopfield, J.J. Neural networks and physical systems with emergent 
collective computational abilities. Proc. Natl. Acad. Sci.USA 79, 2554- 
2558(1982). 

25. Lee, S.H., et al. Efficient generation of midbrain and hindbrain neurons 
from mouse embryonic stem cells. Nat Biotechnol. 18, 675-9 (2000). 

26. Kawasaki, H., et al. Generation of dopaminergic neurons and pigmented 
epithelia from primate ES cells by stromal cell-derived inducing activity. Proc 
Natl Acad Sci U SA.,99, 1580-5 ( 2002). 

27. Westmoreland, J.J. Hancock, C.R. and Condie, B.G. Neuronal 
development of embryonic stem cells: a model of GABAergic neuron 
differentiation. Biochem Biophys Res Commun. 284, 674-80 (2001 ). 

28. Chang, J.C., Brewer, G.J. & Wheeler, B.C. Modulation of neural network 
activity by patterning. Biosensors & Bioelectronics 16 527-533 (2001) 

29. Nakamura, F., Kalb, R.G., & Strittmatter, S.M. Molecular basis of 
semaphorin-mediated axon guidance. J. Neorubiol. 44, 219-229 (2000). 
30Tessier-Lavigne, M. & Goodman, C.S. The molecular biology of axon 
guidance. Science 274, 1 123-1 133 (1996). 

31. Hubel DH, Wiesel TN. Early exploration of the visual cortex. Neuron, 20, 
401-412(1998). 

32. Malgaroli, A. & Tsien, R.W. Glutamate-induced long term potentiation of 
the frequency of miniature synaptic currents in cultured^jjpjQ^a^aJ^eurons. 



Nature 257,134-139 (1992). 




per se e per gll altrl 
Olga Capasso 
(N* d'tecr. 820 B) 



RM20D2 

RIVENDICAZIONI 

1. Dispositivo per la mappatura e il processamento in parallelo di immagini 
comprendente almeno una base "multl electrode array" MEA, su cui e 
cresciuta una coltura omogenea di cellule neuronali, in cui detta base MEA e 
in grado di stimolare e registrare I'attivita elettrica di dette cellule neuronali. 

2. Metodo per il processamento in parallelo di un immagine digitale 
comprendente le seguenti fasi: 

a) mappatura di una immagine digitale (li. 2 (x,y)) (INPUT) con 1 o 2 bit di 
risoluzione (li(x,y)) o l 2 (x,y)) a seconda che si abbia un'immagine 
rispettivamente di 1 o 2 bit) di NxN pixel in impulsi di voltaggio di 2 o 4 livelli di 
intensita applicati ad una matrice NxN di elettrodi integrati su un piattino 
(MEA) su cui sono mantenute in coltura cellule neuronali connesse 
spontaneamante fra loro a formare una rete; 
' b) elaborazione dell'immagine da parte della coltura neuronale con il nucleo o 
kernel di convoluzione: 

hOo.a.t) = A(t) exp((p - p(t)) / 2 a(t) 2 ) (1 ) 

p 2 =x 2 +y 2 

c) registrazione dell'attivita elettrica della coltura neuronale a mezzo di 
registrazione con la matrice di multielettrodi dei segnali elettrici (in voltaggio) 
extracellulari e rivelazione, per ogni singolo elettrodo e in successivi intervalli 
di tempo, dei segnali o spari associati alia generazione dei potenziali di azione 
dei neuroni in coltura. 

3. Metodo secondo la rivendicazione 2 in cui si registra la frequenza di sparo 
FR(x,y,t) (OUTPUT) misurata dall'elettrodo in posizione (x,y) ed in un 
intervallo di tempo centrato in t. 
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4. Metodo secondo la rivendicazione 3 in cui INPUT ed OUTPUT sono legati 
fra loro dall'equazione: 

FR(x,y,t) = li, 2 (x,yr h(p,a,t) (2) 
dove ** indica una convoluzione bidimensionale. 

5. Metodo secondo la rivendicazione 2 in cui, se Pimmagine digitate in INPUT 
6 definita con 8 bit (l 8 (x,y)) f essa viene scomposta alternativamente in 4 o 8 
immagini rispettivamente da 2 o 1bit (l m i) dove m 6 uguale rispettivamente a 2 
o 1 secondo I'equazione: 

8/m 

l8(x l y)=iil mi 2 m(| - 1) (3) 
1 

} 

ed ogni singola immagine l mi viene processata tndipendentemente e poi 
ricomposta in una unica immagine ad 8 bit, in cui Tintero processo 
(scomposizione, filtraggio, ricomposizione) e descritto dairequazione: 

8/m 

Ii lmi 2 m(M) ** h(p,a,t) (4) 
1 

it cui risultato finale e I'immagine l 8 (x,y) processata con una risoluzione ad 8 
bit. 
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